import w2
data_columns = [
'Air Temperature',
'Dew Point Temperature',
'Wind Speed',
'Wind Direction',
'Cloudiness',
'Solar Radiation'
]
ylabels = [
'Air Temperature ($^oC$)',
'Dew Point Temperature ($^oC$)',
'Wind Speed (m/s)',
'Wind Direction (radians)',
'Cloudiness (fraction)',
'Solar Radiation ($W/m^2$)'
]
met_wb1 = w2.read('BerlinMilton2006/2006_Met.npt', 2006, data_columns)
met_wb2 = w2.read('BerlinMilton2006/2006_Met1.npt', 2006, data_columns)
display(met_wb1)
display(met_wb2)
w2.multiplot(met_wb1, ylabels=ylabels, colors=w2.k2)
w2.multiplot(met_wb2, ylabels=ylabels, palette='colorblind')
| Air Temperature | Dew Point Temperature | Wind Speed | Wind Direction | Cloudiness | Solar Radiation | |
|---|---|---|---|---|---|---|
| Date | ||||||
| 2006-01-01 00:00:00 | 0.6 | -0.6 | 3.89 | 4.01 | 10.0 | 0 |
| 2006-01-01 01:00:00 | 0.6 | -0.6 | 3.89 | 4.01 | 10.0 | 0 |
| 2006-01-01 02:00:00 | 1.1 | -0.6 | 3.26 | 4.01 | 10.0 | 0 |
| 2006-01-01 03:00:00 | 1.1 | -0.6 | 3.96 | 4.04 | 10.0 | 0 |
| 2006-01-01 04:00:00 | 1.1 | -0.6 | 3.58 | 4.16 | 10.0 | 0 |
| ... | ... | ... | ... | ... | ... | ... |
| 2006-12-31 20:00:00 | 11.0 | 7.4 | 7.09 | 2.44 | 10.0 | 0 |
| 2006-12-31 21:00:00 | 10.6 | 8.5 | 7.04 | 2.44 | 10.0 | 0 |
| 2006-12-31 22:00:00 | 11.1 | 9.5 | 8.34 | 2.47 | 10.0 | 0 |
| 2006-12-31 23:00:00 | 11.3 | 9.6 | 5.12 | 2.80 | 10.0 | 0 |
| 2007-01-01 00:00:00 | 11.3 | 9.6 | 5.12 | 2.80 | 10.0 | 0 |
8761 rows × 6 columns
| Air Temperature | Dew Point Temperature | Wind Speed | Wind Direction | Cloudiness | Solar Radiation | |
|---|---|---|---|---|---|---|
| Date | ||||||
| 2006-01-01 00:00:00 | 0.6 | -0.6 | 3.89 | 4.01 | 10.0 | 0 |
| 2006-01-01 01:00:00 | 0.6 | -0.6 | 3.89 | 4.01 | 10.0 | 0 |
| 2006-01-01 02:00:00 | 1.1 | -0.6 | 3.26 | 4.01 | 10.0 | 0 |
| 2006-01-01 03:00:00 | 1.1 | -0.6 | 3.96 | 4.04 | 10.0 | 0 |
| 2006-01-01 04:00:00 | 1.1 | -0.6 | 3.58 | 4.16 | 10.0 | 0 |
| ... | ... | ... | ... | ... | ... | ... |
| 2006-12-31 20:00:00 | 11.0 | 7.4 | 7.09 | 2.44 | 10.0 | 0 |
| 2006-12-31 21:00:00 | 10.6 | 8.5 | 7.04 | 2.44 | 10.0 | 0 |
| 2006-12-31 22:00:00 | 11.1 | 9.5 | 8.34 | 2.47 | 10.0 | 0 |
| 2006-12-31 23:00:00 | 11.3 | 9.6 | 5.12 | 2.80 | 10.0 | 0 |
| 2007-01-01 00:00:00 | 11.3 | 9.6 | 5.12 | 2.80 | 10.0 | 0 |
8761 rows × 6 columns
Input files:
data_columns = ['Inflow (QIN)']
qin_br1 = w2.read('BerlinMilton2006/2006_DeerCrk_Qin.npt', 2006, data_columns)
qin_br2 = w2.read('BerlinMilton2006/2006_WillowCrk_Qin.npt', 2006, data_columns)
qin_br3 = w2.read('BerlinMilton2006/2006_IslandCrk_Qin.npt', 2006, data_columns)
qin_br4 = w2.read('BerlinMilton2006/2006_MillCrk_Qin.npt', 2006, data_columns)
ylabel = 'Inflow ($m^3/s$)'
w2.plot(qin_br1, ylabel=ylabel)
w2.plot(qin_br2, ylabel=ylabel)
w2.plot(qin_br3, ylabel=ylabel)
w2.plot(qin_br4, ylabel=ylabel)
qin_br1
| Inflow (QIN) | |
|---|---|
| Date | |
| 2006-01-01 00:00:00 | 1.02 |
| 2006-01-01 01:00:00 | 0.92 |
| 2006-01-01 02:00:00 | 0.91 |
| 2006-01-01 03:00:00 | 1.00 |
| 2006-01-01 04:00:00 | 1.00 |
| ... | ... |
| 2006-12-31 20:00:00 | 1.07 |
| 2006-12-31 21:00:00 | 0.98 |
| 2006-12-31 22:00:00 | 1.07 |
| 2006-12-31 23:00:00 | 1.26 |
| 2007-01-01 00:00:00 | 1.36 |
8761 rows × 1 columns
Two outlets at Berlin Dam
Input files:
data_columns = ['Dam Outflow (QOT): Outlet 1', 'Dam Outflow (QOT): Outlet 2']
qot_br1 = w2.read('BerlinMilton2006/2006_BerlinDam.npt', 2006, data_columns)
qot_br2 = w2.read('BerlinMilton2006/2006_MiltonDam.npt', 2006, data_columns)
ylabel = 'Dam Outflow ($m^3/s$)'
w2.plot(qot_br1, title='2006_BerlinDam.npt', ylabel=ylabel)
w2.plot(qot_br2, title='2006_MiltonDam.npt', ylabel=ylabel)
Input files:
data_columns = ['Tributary Inflow (QTR)']
qtr_tr1 = w2.read('BerlinMilton2006/2006_BeechCrk_Qtr.npt', 2006, data_columns)
qtr_tr2 = w2.read('BerlinMilton2006/2006_Alliance_Qtr.npt', 2006, data_columns)
qtr_tr3 = w2.read('BerlinMilton2006/2006_MahoningR_Qtr.npt', 2006, data_columns)
ylabel = 'Tributary Inflow ($m^3/s$)'
w2.plot(qtr_tr1, title='2006_BeechCrk_Qtr.npt', ylabel=ylabel)
w2.plot(qtr_tr2, title='2006_Alliance_Qtr.npt', ylabel=ylabel)
w2.plot(qtr_tr3, title='2006_MahoningR_Qtr.npt', ylabel=ylabel)
Input files:
data_columns = ['Inflow Temperature (TTR)']
ttr_tr1 = w2.read('BerlinMilton2006/2006_MahoningR_Ttr5.npt', 2006, data_columns)
ttr_tr2 = w2.read('BerlinMilton2006/2006_Alliance_Ttr.npt', 2006, data_columns)
ttr_tr3 = w2.read('BerlinMilton2006/2006_MahoningR_Ttr.npt', 2006, data_columns)
ylabel = 'Inflow Temperature ($^oC$)'
w2.plot(ttr_tr1, title='2006_MahoningR_Ttr5.npt', ylabel=ylabel)
w2.plot(ttr_tr2, title='2006_Alliance_Ttr.npt', ylabel=ylabel)
w2.plot(ttr_tr3, title='2006_MahoningR_Ttr.npt', ylabel=ylabel)
Input files:
infiles = [
'2006_MahoningR_Ctr1.npt',
'2006_AllianceWwtp_Ctr.npt',
'2006_MahoningR_Ctr.npt'
]
ctr1_vars = ['TDS', 'SO4', 'Cl', 'ISS', 'OP', 'NH4', 'NOx', 'Fe',
'LDOM', 'RDOM', 'LPOM', 'RPOM', 'BG', 'DIAT', 'OTH', 'DO']
ctr2_vars = ['TDS', 'ISS', 'OP', 'NH3', 'NOx', 'LDOM', 'RDOM', 'LPOM', 'RPOM', 'DO']
ctr3_vars = ['TDS', 'SO4', 'Cl', 'ISS', 'OP', 'NH4', 'NOx', 'Fe',
'LDOM', 'RDOM', 'LPOM', 'RPOM', 'BG', 'DIAT', 'OTH', 'DO']
ctr1_ylabels = [
'TDS (mg/L)', 'SO4 (mg/L)', 'Cl (mg/L)', 'ISS (mg/L)', 'OP (mg/L)', 'NH4 (mg/L)', 'NOx (mg/L)', 'Fe (mg/L)',
'LDOM (mg/L)', 'RDOM (mg/L)', 'LPOM (mg/L)', 'RPOM (mg/L)', 'BG (mg/L)', 'DIAT (mg/L)', 'OTH (mg/L)', 'DO (mg/L)']
ctr2_ylabels = ['TDS (mg/L)', 'ISS (mg/L)', 'OP (mg/L)', 'NH3 (mg/L)', 'NOx (mg/L)',
'LDOM (mg/L)', 'RDOM (mg/L)', 'LPOM (mg/L)', 'RPOM (mg/L)', 'DO (mg/L)']
ctr3_ylabels = [
'TDS (mg/L)', 'SO4 (mg/L)', 'Cl (mg/L)', 'ISS (mg/L)', 'OP (mg/L)', 'NH4 (mg/L)', 'NOx (mg/L)', 'Fe (mg/L)',
'LDOM (mg/L)', 'RDOM (mg/L)', 'LPOM (mg/L)', 'RPOM (mg/L)', 'BG (mg/L)', 'DIAT (mg/L)', 'OTH (mg/L)', 'DO (mg/L)']
ctr1 = w2.read(f'BerlinMilton2006/2006_MahoningR_Ctr1.npt', 2006, ctr1_vars)
ctr2 = w2.read(f'BerlinMilton2006/2006_AllianceWwtp_Ctr.npt', 2006, ctr2_vars)
ctr3 = w2.read(f'BerlinMilton2006/2006_MahoningR_Ctr.npt', 2006, ctr3_vars)
w2.multiplot(ctr1, ylabels=ctr1_ylabels, figsize=(15, 30), style='.-')
w2.multiplot(ctr2, ylabels=ctr2_ylabels, figsize=(15, 19), style='.-')
w2.multiplot(ctr3, ylabels=ctr3_ylabels, figsize=(15, 30), style='.-')
Input files:
data_columns = ['Distributed Tributary Inflows (QDT)']
qdt_br1 = w2.read('BerlinMilton2006/2006_Bqwb_72smooth.npt', 2006, data_columns)
qdt_br5 = w2.read('BerlinMilton2006/2006_Mqwb_72smooth.npt', 2006, data_columns)
w2.plot(qdt_br1, title='2006_Bqwb_72smooth.npt', xlabel='Date', ylabel=f'Distributed Tributary Inflows ($m^3/s$)')
w2.plot(qdt_br5, title='2006_Mqwb_72smooth.npt', xlabel='Date', ylabel=f'Distributed Tributary Inflows ($m^3/s$)')
Input files:
data_columns = ['Distributed Tributary Inflow Temperature (TDT)']
tdt_br1 = w2.read('BerlinMilton2006/2006_Btdt.npt', 2006, data_columns)
tdt_br5 = w2.read('BerlinMilton2006/2006_Mtdt.npt', 2006, data_columns)
w2.plot(tdt_br1, title='2006_Btdt.npt', ylabel=f'Distributed Tributary Inflow Temperature ($^oC$)')
w2.plot(tdt_br5, title='2006_Mtdt.npt', ylabel=f'Distributed Tributary Inflow Temperature ($^oC$)')
Input files:
BR5: 2006_MahoningR_Ctr2.npt
Input files:
data_columns = ['Precipitation']
pre_br1 = w2.read('BerlinMilton2006/2006_Pre.npt', 2006, data_columns)
pre_br2 = w2.read('BerlinMilton2006/2006_Pre1.npt', 2006, data_columns)
pre_br3 = w2.read('BerlinMilton2006/2006_Pre2.npt', 2006, data_columns)
pre_br4 = w2.read('BerlinMilton2006/2006_Pre3.npt', 2006, data_columns)
pre_br5 = w2.read('BerlinMilton2006/2006_Pre4.npt', 2006, data_columns)
ylabel = 'Precipitation (mm)'
w2.plot(pre_br1, title='2006_Pre.npt', ylabel=ylabel)
w2.plot(pre_br2, title='2006_Pre1.npt', ylabel=ylabel)
w2.plot(pre_br3, title='2006_Pre2.npt', ylabel=ylabel)
w2.plot(pre_br4, title='2006_Pre3.npt', ylabel=ylabel)
w2.plot(pre_br5, title='2006_Pre4.npt', ylabel=ylabel)
data_columns = ['Precipitation Temperature (TPR)']
tpr_br1 = w2.read('BerlinMilton2006/2006_PreT_noneg.npt', 2006, data_columns)
tpr_br2 = w2.read('BerlinMilton2006/2006_PreT_noneg1.npt', 2006, data_columns)
tpr_br3 = w2.read('BerlinMilton2006/2006_PreT_noneg2.npt', 2006, data_columns)
tpr_br4 = w2.read('BerlinMilton2006/2006_PreT_noneg3.npt', 2006, data_columns)
tpr_br5 = w2.read('BerlinMilton2006/2006_PreT_noneg4.npt', 2006, data_columns)
ylabel = 'Precipitation Temperature ($^oC$)'
w2.plot(tpr_br1, title='2006_PreT_noneg.npt', ylabel=ylabel)
w2.plot(tpr_br2, title='2006_PreT_noneg1.npt', ylabel=ylabel)
w2.plot(tpr_br3, title='2006_PreT_noneg2.npt', ylabel=ylabel)
w2.plot(tpr_br4, title='2006_PreT_noneg3.npt', ylabel=ylabel)
w2.plot(tpr_br5, title='2006_PreT_noneg4.npt', ylabel=ylabel)
# Constituents list: Include an extra "constituent" to the extra column that will be added due to a trailing comma.
constituents = ['TDS', 'G1_SO4', 'G2_Cl', 'ISS1', 'PO4', 'NH4', 'NO3', 'FE', 'LDOM', 'RDOM', 'LPOM', 'RPOM', 'BG', 'DIAT', 'OTH', 'DO']
cwo_37 = w2.read('BerlinMilton2006/cwo_37_wdo.csv', 2006, constituents)
cwo_75 = w2.read('BerlinMilton2006/cwo_75_wdo.csv', 2006, constituents)
ylabels = [
'TDS (mg/L)', 'G1_SO4 (mg/L)', 'G2_Cl (mg/L)', 'ISS1 (mg/L)', 'PO4 (mg/L)', 'NH4 (mg/L)', 'NO3 (mg/L)', 'FE (mg/L)',
'LDOM (mg/L)', 'RDOM (mg/L)', 'LPOM (mg/L)', 'RPOM (mg/L)', 'BG (mg/L)', 'DIAT (mg/L)', 'OTH (mg/L)', 'DO (mg/L)']
w2.multiplot(cwo_37, figsize=(15, 30), ylabels=ylabels)
w2.multiplot(cwo_75, figsize=(15, 30), ylabels=ylabels)
# Constituents list: Include an extra "constituent" to the extra column that will be added due to a trailing comma.
constituents = ['TKN', 'TP', 'CHLA', 'TSS']
dwo_37 = w2.read('BerlinMilton2006/dwo_37_wdo.csv', 2006, constituents)
# The following file has times only, no data
# dwo_75 = w2.TimeSeries('BerlinMilton2006/dwo_75_wdo.csv', 2006, constituents)
ylabels = ['TKN (mg/L)', 'TP (mg/L)', 'CHLA (mg/L)', 'TSS (mg/L)']
w2.multiplot(dwo_37, figsize=(15, 15), ylabels=ylabels)
display(dwo_37)
| TKN | TP | CHLA | TSS | |
|---|---|---|---|---|
| Date | ||||
| 2006-01-01 00:00:00 | 0.0000 | 0.0000 | 0.0000 | 0.0000 |
| 2006-01-01 01:00:00 | 1.1375 | 0.0636 | 4.4195 | 7.2167 |
| 2006-01-01 02:00:00 | 1.1373 | 0.0636 | 4.4188 | 7.2096 |
| 2006-01-01 03:00:00 | 1.1370 | 0.0636 | 4.4178 | 7.1997 |
| 2006-01-01 04:00:00 | 1.1369 | 0.0636 | 4.4177 | 7.1983 |
| ... | ... | ... | ... | ... |
| 2006-06-25 17:00:00 | 0.5503 | 0.0479 | 2.0518 | 0.7734 |
| 2006-06-25 18:00:00 | 0.5504 | 0.0479 | 2.0502 | 0.7729 |
| 2006-06-25 19:00:00 | 0.5504 | 0.0479 | 2.0493 | 0.7725 |
| 2006-06-25 20:00:00 | 0.5501 | 0.0479 | 2.0546 | 0.7726 |
| 2006-06-25 21:00:00 | 0.5495 | 0.0479 | 2.0687 | 0.7733 |
4222 rows × 4 columns